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We have used the method of generating functional in imaginary time to derive the current- voltage 
characteristics of a tunnel junction with arbitrary tunneling conductance, connected in series with an 
external impedance and a voltage source. We have shown that via the renormalized charging energy 
and the renormalized environment conductance, our nonperturbative expressions of the total action 
can be mapped onto the corresponding perturbative formulas. This provides a straightforward 
way to go beyond the perturbation theory. For the impedance being a pure resistance, we have 
calculated the conductance for various voltages and temperatures, and the results agree very well 
with experiments. 
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I. INTRODUCTION 

By decreasing the capacitance C of a double tunnel 
junction structure, the role of Coulomb charging energy 
e 2 /2C becomes more important. For example, if C is of 
the order 1 fF and the temperature T is lower than 1 
K, the elementary Coulomb charging energy e 2 /2C dom- 
inates the charge transport through a tunnel junction sys- 
tem, provided that the tunneling resistance Rt is larger 
than the resistance quantum Rx—h/e 2 . The technology 
nowadays can fabricate tunnel junctions with C less than 
1 fF, and hence phenomena related to Coulomb blockade 
(CB) have been extensively studied in recent yearaltj. 

For a tunnel junction system with Rt much larger 
than RK=h/e 2 , the perturbation theory of CB has been 
well established, where one performs low order pertui 
tion expansion in terms of the tunneling HamiltonianEf 
The so-obtained theoretical results are in good agree- 
ment with, experimental observations in multijunction 
structureaTQ. On the other hand, for a single tunnel 
junction, one has to take into account the effect of circuit 
elements which are attached to the single junction in or- 
der to perform measurements. If the circuit impedance is 
much smaller than the quantum resistance Rk, for exam- 
ple, the free space impedance iJ/~377 tt, the CB effect 
can hardly be detected. However, the CB phenomena 
emerge when a large impedance is inserted into pthe. cir- 
cuit in series connection with the single junctioroQtl. In 
a multijunction system, the effect of circuit impedance is 
secondary because there are more than one large junction 
resistance. The influence of circuit elements is a specific 
type of electromagnetic environment (EME) effect. 

In the weak tunneling regime, Devoret et aZuhave used 
the Fermi's golden rule to treat the tunneling Hamil- 
tonian, and have developed a so-called P{E) theory to 
study the Coulomb charging effect of a single-tunnel junc- 
tion surrounded by an EME. Girvin et alu have used 
the Green's function technique to reach the same the- 
ory. They studied a tunnel junction connected in series 
with an external impedance and a voltage source. The 



Hamiltonian of the device consists of four parts 

H = H C + H cl + H t + H cx . (1) 

When a charge Q is added to the junction with junction 
capacitance C, the Coulomb charging energy is 



H c = Q 2 /2C. 



(2) 



The quasi-particle states in the left and the right elec- 
trode are labeled by a for the spin and the transverse 
motion. The longitudinal wave vectors in the right elec- 
trode are represented by k, and in the left electrode by 
q. In standard quasi-particle notations, the Hamiltonian 
H e i for the electrodes has the form 



A: rr 



E 

qa 



£qcr c qa c qa i 



(3) 



Since a is conserved during the tunneling through a junc- 
tion, we use a to specify the tunneling channel. For a 
fixed tunneling channel, let tk qa be the matrix element 
for tunneling from the state qa in the left electrode to 
the state ka in the right electrode. Then, the tunneling 
Hamiltonian is expressed as 



Ht = ^[tkqaCkaCqa exp(-i(p) + H.C.] , 



(4) 



kq<7 



where the operator ip is conjugate to Q via the commuta- 
tion relation [tp, Q]=ie. We should mention that here we 
have neglected the influences of the tunnelinglime, which 
is appropriate for metallic tunnel junctionsocl. Since the 
energies of tunneling electrons are near the Fermi energy, 
tkqa can be well approximated by a constant value t. In 
the presence of a bias voltage V, thg^jxternal impedance 
is simulated by a set of LC circuitsBQ&CrEj as 



Hp* = 



m—l 



2C n 



(eVt- ip - ip m ) 2 
2e 2 L rn 



(5) 



I 



where the operators (p m and q n obey the commutation 
relation [<p m , q n ]=ieSm,n- 

As we mentioned before, in the regime of weak tunnel- 
ing, perturbation theory has been applied to the Hamil- 
tonian Eq. (Q) to, calculate the current as a function of 
the bias voltagea'trta. Beyond the weak tunneling r.e£ime, 
following the pioneering work of Ambegaokar et aZlid, the 
nonperturbative approach of path intcgxal has been used 
to investigate rJthe statistical propertiesE3 and the trans- 
port behavioru in the linear-response regime of tunnel 
junctions (also see Ref. || for the current-bias case). Re- 
cently, using the path integral approach formulated in 
imaginary time to calculate the linear conductances, we 
found that the strong tunneling processes renormalize 
both the charging energy and the external impedance. 
Furthermore, within this imaginary time formalism, we 
have shown that at sufficiently high temperatures, the ac- 
tion of the tunnel junction can be approximated by the 
action of the corresponding resistance, and so an explicit 
expression of the high temperature quantum conductance 
can be derivedEI. Along the line of Feynman's formalism 
for polaron, the Coulomb charging effectias been inves- 
tigated with the Ohmic approximationEia, the nonequi- 
librium_.Green's function methodli-3, and the mean-field 
theoryEj. The mean-field theory also demonstrates the 
normalization of the charging energy and the external 
impedance, and the so-calculated zero-bias conductance 
agrees with experiment as long as fc^T is not much less 
than e 2 /2C and Rt is not much smaller than Rk- 

Nonperturbative approach has so far produced con- 
crete results only for the situation of zero-bias. In this 
paper we will consider finite bias and so investigate the 
full current-voltage (I-V) characteristics of a voltage- 
biased tunnel junction connected in series with an exter- 
nal impedance. The derivation of the path integral repre- 



sentation for the generating functional in imaginary time 
will be outlined in Sec. II, leaving the complete mathe- 
matical manipulations in the Appendix. With the gen- 
crating functional we obtain the tunneling current as a 
function of the bias voltage. Specific cases suitable for 
comparing with experiments will be studied in details in 
Sec. Ill, showing very good agreement with experiments. 
Finally, in Sec. IV we will prove that via the renormalized 
charging energy and the renormalized environment con- 
ductance, our nonperturbative expressions of the total 
action can be mapped onto the well-established pertur- 
bative results. This provides a straightforward way to go 
beyond the perturbation theory. 

II. GENERAL FORMULATION OF THE I-V 
CHARACTERISTICS 

To study the physical system described by Eqs. (|l|)- 
(|5|), one needs to deal with the electron tunneling and 
the effect of external impedance. While the EME has 
been treated nonperturbatively by many authors, in this 
paper we will analyze both the tunneling part and the 
EME part with a nonperturbative approach. To achieve 
this goal, we will start from the generating functional for 
a tunneling junction connected in series with an external 
impedance and a voltage source. 

The effect of the bias voltage V can be incorporated 
into the Hamiltonian not explicitly depending on the time 
variable via a time-dependent unitary transformation de- 
fined as 

U(t) = ex V [ieVt(J2cl a c ka )}. 

ka 

It transforms the Hamiltonian in Eq. (fy) to 



H = Q 2 /2C + ^2(s krT + eV)c\ a c kcr + e q<y c\ a c qa + [t kqa c\ a c q<J exp (i<f>) + H.c] 

k,a f?,c kqa 



E 



2C n 2e 2 -L n 



(6) 



where Q=Q-CV reflects the quantum fluctuations of the 
charge on the tunnel junction, and <fi=eVt-tp serves as the 
phase of the external impedance. These two new canoni- 
cal variables satisfy the commutation relation [</>, Q]=ie. 
The corresponding current operator has the form 



The generating functional in imaginary time is defined 



as 



Z v [rj\ = tr i T T exp i - J dr [H - / r 7?(r)] \\ , (8) 



I T = -»e y^fcggcj^Cgg exp (i<f>) - H.c] 



kqa 



(7) where T T is the time-ordering operator in imaginary time. 
Substituting 
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U - Itv{t) = Q 2 /2C + J2( £ ka + eV)cl a c ka + ^ e qa c\ a c qa 



^2{tkqa [1 + iev(r)]cl a c qa exp {i<j>) + H.c.} + 

kqa n— 1 



2C n 2e^L n 



(9) 



into Eq. (||), we can rewrite Zv[rj\ as a series expansion 
in terms of the tunneling part of the quantity Ti-Ir'ni' 1 ')- 
Then, (3 is divided into segments of length e — ► 0, and 
eigenstates of </> and Q are inserted between the elements 



of the Trotter product. Next, the charge numbers are 
summed over. Such mathematical manipulations are per- 
formed in the Appendix, leading to the expression 



Z v [rj\ =AA"^(-l) ; f dn f 'dn-i--- f 2 d Tl JT f D^ m f d^f[ f dfr V] ••■ V 

in JO JO JO „ , J J • , J , , , , 



M 



exp{- ^ e{tf/4E c + C m ^j2e 2 + (<p m - faf /2e 2 L m ]} 

xi^i6 • • - + *eCi»7(7Tt)] ■••[! + ^(n)] n^ 6 "^'^.^.'.^^) ' ' ' K^^M"' 



(10) 



where 



/V = ^(E/ccr + e^)c+ r C fecr + ^ EgaC^Cqa (11) 



and 



(12) 



Here we have introduced the notations cf =ci" for £=1, 

ka ka = ' 



and c| =Cfc CT for £=-1. The same notations apply to (A 



qa- 



The trace over the electronic states in Eq. (|io| ) will 
be evaluated by first introducing the quantum statistic 
average of the products of creation and annihilation oper- 
ators, and then using the Wick's theorem. Here we have 
made use of the fact that the number of tunneling chan- 
nels in metallic junctions, N is very large. Therefore, 
terms of order 1/N are neglected. Again, these algebraic 
works arc done in the Appendix, and the path integral 
representation of the generating functional is finally de- 
rived as 



ZVM = II / D fm / DtfKTh ^/ i E c+ C m rtj2eZ + (^ m r / 2e*L m] 

00 i ( fP fp" 1 r 



r=0 



JO 



which is Eq. (15) in the Appendix. 

In the above generating functional, the path integrals 
over the EME mades are Gaussian, and thus can be car- 
ried out exactlyt^l to give 



Z V [ V } = J Dcbe- 



-Sex[0]-S t [0,7,] 



(14) 



In the total action 

S&rjj^SM+S^+St&v], (15) 
the first term on the right hand of the above equation 



SM = 



4£, 



drcj) 2 



c Jo 



(16) 



is resulted from the Coulomb charging energy, the second 
term 

SexM = o/ dT / dT'a C x(T-T')[0(r)-0(r')] 2 (17) 
* Jo Jo 

is due to the external impedance, and the third term 



rP rP 
/ dT I dr'a^T — T 1 ) 
Jo Jo 
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x exp[(eVY + i0(t)) - (eVr' + «<^(r'))] 

[l + ie*7(r)][l-ie»7(T / )] (18) 

describes the contribution of tunneling processes. In 
Eq. (0), a cx (uji)=R K \uji\/4irZ cx (-i\uji\) and Z cx (cj) is 
the Fourier coefficient of the external impedance, where 
u>i=2irl/(3 is the Matsubara frequency. For a purely re- 
sistive impedance i? x, this expression is simplified to 
a cx (aji)=-a cx \u!i\/4Tr with a cx =R K /R cx - The S t [(f>,r)} 
in Eq. (|lj) contains both the bias voltage and the driv- 
ing source of the generating functional. Here at(tOi)=- 
a t \u!i\/4:ir with a t =RK / 'Rt- It is important to notice 
that for V=0 and r](T)=r)(r')=0, Eq. Q) reduces to the 
familiar form of the partition function of an unbiased 
deviceO. 

Before analyzing the current, we would like to clar- 
ify what is indicated by the equations of the generat- 
ing functional and the related actions. In derivation of 
the generating functional Eq. (|lj), the technique that we 
used singles out the ground state, projecting out all ex- 
cited states. At sufficiently low temperatures, our calcu- 
lation should reproduce expectation values in the ground 
state. The ground state can be obtained from the non- 
equilibrium state at finite voltages by the transfer of a 
macroscopic number of electrons from one electrode to 
another. The exponential dependence on the voltage fa- 
vors paths in which a large number of electrons move 
from one electrode to the other at the beginning of the 
path, as implicit in the action given by Eq. JT§|). 

The dc current I(V) is readily derived as the first order 
functional derivative of Zy[r)] with respect to r/, 



I(V) = 2e J daati^Zy 1 J Dcfx 



-s[<p] 



- exp (-eVa) exp - 0(0)])} , (19) 

where Zy is the path integral with the action 
S[4>]=S Q [4>]+S e x[4>]+S t [(t>,r)=0]. The current auto- 
correlation function and high order correlation functions 
can be obtained as well with the corresponding high order 
functional derivatives of the generating functional Zy [r/] . 

Since in standard four-probe experiments, the tunnel- 
ing current is measured as a function of the average volt- 
age drop Vf across the tunnel junction, instead of the 
bias voltage V, we should rewrite the above equation in 
a suitable form. This can be done by introducing the 
phase fluctuation of the external impedance, 8(t)=<f>(t)- 
e(V — Vt)t, into the path integral expression for the tun- 
neling current. Then, Eq. (091) becomes 



/3 



I = 2e / daa t (cr)Z- L / D9e 



x - {exp (eV t a) exp(i[9(a) - 6(0)}) 
- exp (-eV t a) exp (-i[6(a) - 9(0)})} 



where Z is the partition function. The total action 
S\9}=S c \6}+S t \0}+S cx \8} in the above equation is read- 
ily obtained from Eq. ( |15| ) by setting 77=0 and V=0, but 
with 8 instead of tfi as the variable. 

Before going further, we would like to point out th 
the above general expression reduces to the P(E) theor 
at the weak tunneling limit. To show this, we only need 
to neglect the contribution of the tunnel junction to the 
total action, and then to perform an analytical contin- 
uation from imaginary time to real time. Therefore, for 
tunnel junctions with tunnel resistances much larger than 
the quantum resistance, one can either use the P(E) the- 
ory formulated in real time, or use our Eq. ( p0[ ) repre- 
sented in imaginary time. 

In experiments, it is convenient to use Cr<-£k«s near 
the tunnel junction as well-controllable EMEtaO. Con- 
sequently, from now on we will focus our attention on the 
case of purely resistive EME. Via a series expansion in 
(3eV, we will derive the required linear and nonlinear re- 
sponse functions. We must point out that our results are 
not accurate if j3eV is large. Hence, when we compare 
our theory with experiments in the next section, we have 
checked that the theoretical results are valid under the 
conditions the experiments were performed. 



III. RESULTS COMPARED WITH 
EXPERIMENTS 

In this section, we will derive from Eq. ( |2(i| ) some ex- 
plicit results suitable for comparing with experiments!!-! 
For this purpose we will follow the approach of equivalent 
circuits, where the total action of the EME including the 
tunnel junction is replaced by the actions of an effective 
resistance and an effective capacitance. We will at first 
calculate the effective circuit parameters, and then use 
them to present the full I-V characteristics of the tunnel 
junction. 

We need to evaluate the partition function Z in 
Eq. (20). Let us first write down the well-established 



partition function in the weak tunneling limit 



(oh-i 



where 



^=n(^ (o) ) 

2=1 



(0) _ "cx^l 

1 ~ 2E r 2ir 



(21) 



(22) 



(20) 



are eigenvalues with respect to the eigenfunctions of the 
action 

00 

S o [0] = S c [6} + S cx [9} = J2 + Of) . 

1=1 

Beyond the weak tunneling regime, the partition func- 
tion, although more complicated, can still be derived in 
the same manner. To the first order in at, we have 



4 



z = fii + ^W^V'-\ (23) 



n=l P A i 

with the function /(r) is defined as 

_ /£>gexp{-5o[fl(r)]±z[fl(T)-fl(0)]} 
AJ /^exp{-5 [fl(r)]} ' 

At sufficiently high temperatures, with a Taylor series 
expansion of exp [f(r — t')], the partition function Z can 
be expressed in the same form as Z in Eq. but with 
the eigenvalues 



2£* 



2tt 



(25) 



Here the effective dimensionless conductance a* = a cx 
+ at+0(f3E c ) contains the contributions from both the 
external resistance and the tunnel junction. The effective 
charging energy E*, which is determined from 



1 



f 



I3E* I3E C 



7.2a t f3E c 
8tt 4 



0[(j3E c 



(26) 



depends not only on the bare charging energy E c , but 
also on the junction conductance at- The effect of the 
total EME is then embedded in the effective capaci- 
tance e 2 /2E* and the effective environmental resistance 
Rk/cx* . Consequently, from Eq. (20) the,transport prop- 
erties can be investigated analytically!!, and the total 
conductance 



G{V t ) = G(0) + SG(V t ) 
has a voltage-independent term 



(27) 



G(0) 



f 



{1 



R T 3 



i 

15 



7.2a* 
4vr 4 



*\2 



(28) 



and a low-voltage correction 

+0[((3eV t ) 4 } . 



0[{f3E*f]} 



(29) 



The dots in Fig. 1 are the measured normalized con- 
ductance taken from the Fig. 2a of Ref. ( |l7j ). Using the 
sample parameter values at=1.04 and a e x=8.09, as well 
as the experimental temperature T=4.2 K, our calcu- 
lated normalized conductance as a function of the voltage 
is plotted in Fig. 1 as the solid curve. Our analytical re- 
sult is in very good agreement with the measurement not 
only for the zero- voltage conductance, but also for the en- 
tire curve. We notice that the experimental data exhibits 
conductance-step structures with sharp changes at volt- 
ages V±i~±0.2mV and V±2— ±0.4 mU, a phenomenon 
which was not discussed in the original experimental pa- 
per Ref. ( [l?]). Our conjecture is that such conductance 



steps are originated from the resonances of the EME 
modes. Using a single LC mode as the EME,and in the 
weak tunneling regime, Ingold and Nazarovn have pre- 
dicted the conductance step structures, which are sharp 
at low temperature but are smeared out at higher tem- 
peratures. If we estimate the inductance L from the step 
width \eV±i\—h/\/LC with the capacitance in the range 
of fF, we found L in the range of pH-. which is in agree- 
ment with the experimental valueaHI. These resonant 
EME modes are not included in the purely Ohmic EME, 
but can be investigated if we use the frequency-dependent 
external impedances to model the EME. 

The temperature dependence of the zero-voltage con- 
ductance has been measured thoroughly in Ref. ( p~7| ) . 
For the standard four-terminal setup, if the temperature 
is not too low, the inverse of the conductance dip at zero 
voltage is linear in temperature 

(AG/G T ) -1 = G T /[G T - G(0)] = 3k B T/E c + S . (30) 



Our theory gives the offset 

5 = 0.6 + 0.167(a ox + a t ). 
while the P(E) theory predicts 



0.6 + 0.167a o 



(31) 



(32) 



Two sets of experimental data, taken from the Figs. 3a 
and 3b of Ref. ( [l?]), are plotted in Fig. 2 as dots, to- 
gether with the best fitted dotted lines. The sample pa- 
rameter values are a t =5.86, a cx =20.32 and C =1.99 fF 
for the upper set of data, and at=3.02, a GX =1.50 and 
C=0.92 fF for the lower set of data. Using these param- 
eter values, the analytical results of our theory are plot- 
ted in Fig. 2 as solid lines, and those of the P(E) theory 
are shown in dashed lines. While our theory agrees very 
well with the measurements, the deviation of the P(E) 
theory from the experiments increases when the tunnel- 
ing conductance gets larger. We have_also performed the 
self-consistent numerical calculationsEj formulated in real 
time, and the results are indistinguishable from our ana- 
lytical solution given by Eqs. (|0|) and (|3l|). 



IV. ANALYTICAL SOLUTIONS BASED ON 
RENORMALIZED PARAMETERS 

In the previous section we have shown that in the 
regime of not too low temperatures, our strong tunnel- 
ing formulas can be mapped onto the weak tunneling 
formulas by renormalizing the charging energy and the 
external impedance. Here we will analyze the situation 
of low temperatures. With large (3, the function /(r) in 
Eq. (p3) can be written in a simple analytical form 



f(r) 



2 



7 + In ■ 



(33) 
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where 7 is the Euler constant. With this expression we 
perform a similar mathematical manipulation as what we 
have done for the case of high temperature. For a ex ^>l, 
the so-obtained eigenvalues have again the same form 
as given by Eq. (J2q), but with the effective dimension- 
less conductance a*=a c ^+0(l/(3E c ), and the effective 
charging energy 

47iT(2/a cx ) 

In the above equation T(x) is the Gamma function. 
At low temperatures the renormalization is expected to 
be weak, because the tunneling processes suffer strong 
Coulomb blockade. Using Eqs. @ and @, the I-V 
curves are derived analytically as 

2/Qcx 



v t 



I 



R T T(2 + 2/a e 



7rexp(-7)|eU f | 



(35) 



Now we have derived the effective conductance a* and 
the effective charging energy E* analytically for both low 
and high temperatures. Thus, we can map the nonper- 
turbative expressions of the total action onto the per- 
turbative results by replacing the bare parameters a cx 
and E c with the effective parameters a* and E*, re- 
spectively. Consequently, we can calculate the tunnel- 
ing current and the current-current correlation function 
beyond the perturbation theory. A very important fea- 
ture of our renormalization theory is the appearance of 
the junction conductance in the effective charging en- 
ergy E* as well as in the effective environmental conduc- 
tance a*. Consequently, it is inappropriate to approxi- 
mate the tunnel junction simply by an Ohmic element, 
although it is- a good approximation at sufficiently high 



temperatureaJ'E^E 
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APPENDIX: 

In this Appendix we will derive Eq. (|l3|) for the path 
integral representation of the generating functional of a 
voltage-biased tunnel junction connected in series with 
an external impedance. Since the cases we are interested 
in are not restricted to weak tunneling, we will use the 
nonperturbative resummation technique. 

^From Eqs. (|^) and (^), we can expand the generating 
functional in a Taylor series of the tunneling processes 

00 rfi rri />T2 

ZvM=J2(- i y dTl *n-i~- d n 

l=0 Jo Jo Jo 



tv{e X p[-pn ]n t (n)---n t {T l )} 

where 

H = Q 2 /2C + J2( £ k° + eV)c\ a c k<y 



(Al) 



2C T , 



(0 - 'fin) 



2 1 



2e 2 L r: 



, (A2) 



and 



U t = ^{ifcgtxl 1 + ier l( T )}cia c <i<? ex P + H - c -i • ( A3 ) 

kqa 

Since the EME is modeled by a set of harmonic oscilla- 
tors, the trace over the EME modes can be expressed ex- 
plicitly as path integrals. By dividing j3 into M segments 
of length £=/3/ (M+l) — > and then inserting eigenstates 
of 4> and Q between the elements of the Trotter product, 
we have 



ZvM = ^'£(-1)' J Q dn jf ' dr ; _! ■ • • jf 3 dn fj J Dtp m J dfa J dfa ■ ■ ■ J d<f> M 
J2 ■■■ J2 ^p{~e[E c n 2 + C m up 2 j2e 2 + (<p m - fa) 2 /2e 2 L m }} 

k 1 ,q 1 ,(T 1 ,C 1 k l ,q l ,a l ,C l n 1 n 2 n M+1 

x exp{mi(0 o - fa)} e^p{~e[E c nl + C m ip 2 Jl /2e 2 + (cp m - fa) 2 /2e 2 L m }} exp{in 2 (fa - fa)} ■ ■ ■ 
x ex.p{-e[E c ti 2 M+1 + C m (p 2 n /2e 2 + {ip m - cj> M+ i) 2 /2e 2 L m ]} enp{in M +i{fa\i - fa)} 

xt ! £i6 •••&[! + ieZivin)} •••[! + ie&jjfa)] t^e^/i^^ ( n ) ■ ■ • /^.^(n), (4) 
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where tiv is the trace over all electron states in the chan- 
nel a, 

K = ^2(£ka + eV)c\ a C k<7 + E £ q<? c ta c q<? ( 5 ) 



and 



h M (r) = e^clc^e-^e^ . (6) 



To simplify our mathematical expressions, here we have 
introduced the notations 0^=0^ for £=1, and S k<y =Ck a 
for £=-1. The same notations apply to c q<7 . 

Using the Poisson's resummation formula to sum over 
m, we obtain 



rii — — oo 
oo 

x e 

p— — oo 



7rp4>i/E c -ir 2 p 2 /eE c -e4> 2 /4E C 



(7) 



where (f>i—((f>i-(j)i^i)/e as e — > 0. In the limit sE c — > 0, 
only the p=0 term is relevant and all other terms are 
exponentially small. Then we have 



urn y 



-eE c n i —iriie4>i _ 



n,= — oc 



-ej>1/iE c 



(8) 



After the summations over m have been performed, the 
generating functional is simplified to the form 



Z v [n] = Af"J^(-l) l f dn J ' dn-! ■ ■ ■ J 2 d Tl f[ J D^ m J dfofl J dfa E ••• E 

1=0 m=l " i=l k 1 ,q 1 ,a 1 ,$ 1 k l ,q l ,a l ,£ l 



A I 



exp{- 4ti/^E c + C m <p 2 m /2e 2 + {ip m - /2e 2 L m ]}t l ^ 2 ■ ■ ■ 6 [1 + ie^ivin)} ■■■[1 + ie^v(n) 
»=i 

x[Jtr ff e-^ft liig|i(r|i{i (7j) • ■■/ifc 1 ,, 1 , < r 1 , £l (Ti) . 



(9) 



Now we need to evaluate the trace over the electronic 
states in the above equation. This can be done by divid- 
ing the generating functional with a constant term 



J3tr CT exp[-/3(^ 



which can be absorbed in an irrelevant prefactor TV. In 
this way we have 



ZvH = JVX(-i)' f d n P dn-! ■■■ p dn I ' D<f>f[ J Dy m £ • • • E 

(=0 m=l k 1 ,q 1 ,a 1 ,£ 1 k l ,q l ,a l ,£ 

^16 - 6[i + ietnin)] - [i + i^Me" ^/«.^«./^ + (^-*(r))V^^] 



(10) 



where the symbol (O)o denotes the quantum statistic av- 
erage over free quasi-particles, 



<0>o = 



n.tr 



O 



• (11) 



Since only the combinations (cf ct (t)c. ^(t'))o and 

(c~i(r)c| . CT (r'))o are nonzero, we use the Wick's theo- 
rem to obtain 



E - E 

k 1 ,q 1 ,a 1 ,^ 1 k l ,q l .a l .^ l 



t l ^2 ■ ■ ■ ^'^ ( ^ +eV ^ ■ ■ ■ e «i [x + j^^^)] . . . [i + ie^n) 



x ( c v, fa)^, (n) ■ ■ ■ (n)4; CTi (n)) 



„«1 



7 



1/2 



pairs p=l £ 



xe 1 1 L e 2 



(12) 



where 

a t (r)^ 2 ^ (4 2CT (r) C ^)o(c-i(r)4 i(r )o, (13) 

k 1 ,q 1 ,k 2 ,q 2 ,<7 

and the time-variables r„ and t„ are taken from 
the values {ri, t%, ■ ■ ■ t;}. We notice that the contri- 
butions of the terms with four or more channel in- 
dices equal are of order 1/N. In metallic tunnel junc- 
tions, N is very large and so such terms can be ne- 
glected. Then for nonzero r, the kernel function a t (r) 



and because the sum over r pairs gives a factor (2r — 1)!!, 
we arrive at the path integral representation of the gen- 



can be expressed as a t (r)=a t /4/3 2 sin 2 (ttt / 0) , where 
a t =Rx I 'R,T=4:TT 2 \t\ 2 pip r N with pi (or p r ) being the den- 
sity of states in the left (or right) electrode. For Matsub- 
ara frequencies uj u =2iti>/ (3, the relevant Fourier compo- 
nents reduce to the simple form ott(u>v)=— at|w„|/47r. It 
is obvious that the nonzero contributions to the sum over 
pairs in Eq. ( p"2| ) come from those terms with Yl^-y £i=0 
and even I. We let l=2r be such even integers, and then 
because 



(14) 



erating functional 



r 

II Yl %v-?p 2 a *( T Pi - r P 2 )[ 1 + ie €pMT Pl )][l + *ef p jyfoj] 

= T n a*(r Pl - r Ps )[l + ier,(T Pl )][l - ^^r^le 1 ^^^^^ 1 -^^'^^ 1 , 



m=l J r=0 T ' 

x J jf rfr | <ir'a t (T - /JeK^^W)-!^'^^'))!!! + ie ^ r )][i _ j e?) ( T ')] j . ( 15 ) 



This is Eq. dT|) in Sec. II. 
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FIG. 1. Normalized conductance of the tunnel junction as a function of the junction voltage for a t = 1.04, a cx =8.09 and 
PE C =0.18 at T=4.2 K. The solid curve is from our theory, and the dots are experimental data taken from Fig. 2a of Ref. (17). 
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FIG. 2. Inverse of the normalized conductance dip at zero voltage as a function of the temperature for a t =5.86, a cx =20.32 
and C=1.99 fF (upper curves); and a<=3.02, a cx =1.50 and C=0.92 fF (lower curves). The dashed lines are calculated from 
the P(E) theory, while the solid lines are calculated from our present theory. The dots are the corresponding experimental data 
taken from Figs. 3a and 3b of Ref. (17), together with the best fitted dotted lines. 
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